Glass transition in the quenched and annealed version of the frustrated lattice gas 

model 
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In this paper we study the 3d frustrated lattice gas model 
in the annealed version, where the disorder is allowed to evolve 
in time with a suitable kinetic constraint. Although the model 
does not exhibit any thermodynamic transition it shows a 
diverging peak at some characteristic time in the dynami- 
cal non-linear susceptibility, similar to the results on the p- 
spin model in mean field and Lennard- Jones mixture recently 
found by Donati et al. [cond-mat/9905433]. Comparing 
these results to those obtained in the model with quenched 
interactions, we conclude that the critical behavior of the dy- 
namical susceptibility is reminiscent of the thermodynamic 
transition present in the quenched model, and signaled by 
the divergence of the static non-linear susceptibility, suggest- 
ing therefore a similar mechanism also in supercooled glass- 
forming liquids. 



I. INTRODUCTION 

The study of glass forming systems and spin glasses 
has shown that these systems present a similar complex 
dynamical behavior. In both cases the relaxation time 
increases dramatically when the temperature is lowered; 
furthermore at temperatures lower than some tempera- 
ture T* , the relaxation functions are well fitted at long 
times by a stretched exponential function: 



/(*) = Jo exp 



t 



0</3 < 1. 



(1) 



This similarity was further stressed by the observation 
that the dynamical equations of a class of mean field 
spin glass models, called p-spin glasses [1], are precisely 
equal to the mode coupling equations for supercooled liq- 
uids. The p-spin glasses are a generalization of the spin 
glass model, where spins interact via three or more body 
interactions. 

Despite these analogies the connection in finite dimen- 
sion between glass forming systems and spin glasses is not 
completely clear. As far as the static properties are con- 
cerned, spin glasses undergo a thermodynamic transition 
at a well defined temperature Tsg, where the non-linear 
susceptibility diverges (a similar behavior is found also in 
p-spin in finite dimension [2]). The class of systems that 
show a transition of this kind contains systems with very 
different microscopic structures, but with two essential 



common characteristics, namely the presence of compet- 
itive interactions (frustration) and of quenched disorder. 
On the other hand, glass formers are a class of systems 
where disorder is not originated by some fixed external 
variables, but is "self-generated" by the positions and 
orientation of particles. Differently from spin glasses, in 
glass forming liquids there is no sharp thermodynamic 
transition, characterized by the divergence of a thermo- 
dynamical quantity analogous to the non-linear suscep- 
tibility in spin glasses. However Donati et al. [3] have 
recently introduced a time dependent non-linear suscep- 
tibility both for spin models and for structural glasses: 
they have shown that in the mean field spherical model 
(where the mode coupling equations are exact), and in 
a Lennard- Jones mixture studied by molecular dynam- 
ics simulation [4], the dynamical susceptibility exhibits 
a maximum at some characteristic time, and that this 
maximum diverges as the dynamical temperature Trj is 
approached from above. One may wonder whether the 
presence of this maximum is related somehow to the pres- 
ence of "quasi" -quenched disorder self-generated in the 
liquid, resembling the divergence of the non-linear sus- 
ceptibility present in spin systems with quenched disor- 
der. To shed some light on this problem in this paper 
we want to compare the properties of the frustrated lat- 
tice gas model, which has been recently introduced in 
the context of the glass transition [5], in two cases: when 
the disorder is quenched and when the disorder is self- 
generated. 

In the quenched case the model is a spin glass diluted 
with lattice gas variable, that being constituted by dif- 
fusing particles is suited to study quantities such as the 
diffusion coefficient or the density autocorrelation func- 
tions, that are usually important in the study of liquids. 
This model exhibits in mean field [6] properties closely 
related to those of p-spin models. In 3d, at low enough 
temperature, numerical simulations [7] show a behavior 
of the diffusion coefficient very similar to that experimen- 
tally observed in fragile liquids [8]. Moreover the model 
presents a continuous static transition where the fluctu- 
ations of the order parameter, which coincide with the 
non-linear susceptibility, diverge. This property is ab- 
sent in the infinite dimension p-spin models [3] and in 
supercooled liquids studied by molecular dynamics sim- 
ulations, as Lennard- Jones binary mixtures [4]. 

Here we show the results obtained by studying the 
3d frustrated lattice gas model in the annealed version, 
where the interactions are allowed to evolve with a kinetic 



constraint [10]. We find that the dynamical behavior 
fits quite well the behavior predicted by the mode cou- 
pling theory. It is also easy to show that the model does 
not present any thermodynamic transition, consequently 
there is no divergence in the non-linear susceptibility. 
However the dynamical non-linear susceptibility exhibits 
a maximum with a behavior similar to that found in the 
mean field p-spin model and in molecular dynamics sim- 
ulation of Lennard- Jones mixture. We compare this be- 
havior with the corresponding quantities calculated in 
the quenched version and conclude that the thermody- 
namic transition, present in the quenched model and sig- 
naled by the divergence of the static susceptibility, man- 
ifests itself in the annealed model in the critical behavior 
of the dynamical susceptibility. This behavior seems a 
consequence of the fact at short enough time the interac- 
tions can be considered as quenched variables. Since the 
annealed model shows a behavior reminiscent of super- 
cooled glass forming liquids, these results suggest that 
also in glass forming liquids the behavior of time depen- 
dent non-linear susceptibility can be due to the presence 
of slow degree of freedom which acts for short time as 
quenched variables. Although the annealed lattice gas 
model does not show any thermodynamic transition, we 
cannot exclude that this is due to the absence of signifi- 
cant interactions: typically in a real system the dynam- 
ical constraint corresponds to some kind of interactions; 
in our case the static of the model is instead described 
by a trivial Hamiltonian, while the complex dynamics is 
due to the kinetic constraint. 

In Sect. II we recall briefly the main results of the mode 
coupling theory for supercooled liquids; in Sect. Ill we 
introduce the dynamical non-linear susceptibility as de- 
fined by Donati et al. Finally in Sect. IV we present the 
frustrated lattice gas model and the dynamical behavior 
observed by numerical simulations in the quenched (Sect. 
IV A) and in the annealed version (Sect. IV B). 



nent parameter 0.5 < A < 1, that is a constant depending 
only on the system, and the separation parameter a, that 
is proportional to x — x c , where x is the external control 
parameter (density or temperature) and x c is the critical 
value (c is chosen positive in the glass phase). Via the 
transcendental equation 



r 2 (l-a) _ T 2 (l + 6) 
r(l - 2a) ~ r(l + 26) 



A. 



(2) 



the exponent A determines two exponents, < a < 0.5 
and < b < 1, that rule the relaxation of the system 
near the critical point. 

The MCT predicts that in the j3— regime, near the dy- 
namical transition, the correlators can be written as 

^q(t) = fq + h q c a g±{t/t a ), for t < t < r Q (3) 

where c a = y|a| an d ± refers respectively to the glass 
and liquid phase. The exponent a fixes the short time 
behavior, g±(t/t a ) — (tjt a )~ a for to *C t <C t a , while 
for („ « ( « r a one has a constant in the glass phase 
g+(t/to) = (1 — A) 1 / 2 , and the so called von Schwcidlcr 
law in the liquid phase g-(t/t a ) = —B{t/t a ) h . Here two 
time scales, diverging as the critical point is approached 
from above, appear: 



\-S 



t a = t \<T\ 
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(4) 
(5) 



and to is a microscopic time characteristic of molecules 
motion. 

Furthermore, in the liquid phase, the theory predicts 
the following scaling law for the a— relaxation (t ^> r a ): 



*«(*) 



Mi 



(6) 



II. THE MODE COUPLING THEORY 

In order to compare the dynamical behavior of the frus- 
trated lattice gas with the predictions of the mode cou- 
pling theory (MCT) for supercooled liquids [9] we recall 
briefly the results of this theory. The equations of mo- 
tion of the normalized spatial Fourier transform of the 
density autocorrelation functions, & q (t), are evaluated 
under suitable approximations, and a dynamical transi- 
tion, considered as an idealization of the glass transition, 
is observed: at high temperature the solutions & q (t) van- 
ish at long time (liquid phase); at temperatures below a 
certain critical value Tmct solutions with a nonzero long 
time limit f q (called Debye- Waller factor) appear (glass 
phase). This transition is due to the non linearity of 
the equations, and no thermodynamic phase transition 
is present. 

Let us introduce some important quantities: the expo- 



where the master curve $ g (i/r Q ) is well fitted by a 
stretched exponential of the form (1), with < j3 < 1 
depending on the particular correlator, but not on the 
temperature. This functional form usually fits the ex- 
periments as well. 

In order to test the predictions of the MCT in the 
(3— regime, Gleim and Kob [11] have introduced the fol- 
lowing quantity: 



R q (t) 



<M*)-<M*0 



(7) 



where t' and t" arc arbitrary times in the /3-rcgimc 
(f ^ t"). ^From Eq. (3) we can see that, if <& q (t) is 
in agreement with the leading-order prediction of the 
theory, then R q (t) must be independent on q in the 
0— regime. 



III. THE DYNAMICAL NON-LINEAR 
SUSCEPTIBILITY 

Donati et al. [3] have recently defined a dynamical non- 
linear susceptibility both for spin models and for struc- 
tural glasses. They have shown that in the mean field 
p-spin spherical models (where the mode coupling equa- 
tions are exact) there is a characteristic time where the 
dynamical susceptibility has a maximum, and that this 
maximum diverges as the dynamical temperature To is 
approached from above. 

The Hamiltonian of the p-spin model is 



H = 



E 



'A'i 






(8) 



«i<-"<i p 



where p > 3, the couplings Ji t ...i are Gaussian with zero 
mean and variance 1/N p ~ l , and the spins are real vari- 
ables, with the global constraint X)i=i $f = N, where N 
is the number of spins. The dynamical non-linear sus- 
ceptibility x(t) is defined by 



X (t) = PN((q(t) 2 ) (q(t)) 2 ), 



(9) 



where q{t) = jj X)i Si(t')Si(t' + t) is the overlap between 
the states at times t' and t' + t. Solving the equation 
of motion for x{t) at temperature higher than To, they 
find that \(t) displays a maximum as a function of time, 
x{t*), which is shifted to larger times t* as T approaches 
To from above and increases as a power law x(t*) ex 
(T-T D )- a . 

A similar behavior has also been found in molecular 
dynamics simulation performed for a Lennard- Jones mix- 
ture. 



IV. THE FRUSTRATED LATTICE GAS MODEL 



properties recalling those of a "frustrated" liquid. In fact 
the first term of Hamiltonian (10) implies that two near- 
est neighbor sites can be freely occupied only if their spin 
variables satisfy the interaction, that is if eijSiSj = 1, 
otherwise they feel a strong repulsion. 

To make the connection with a liquid, we note that the 
internal degree of freedom Si may represent for exam- 
ple internal orientation of a particle with non symmetric 
shape. Two particles can be nearest neighbors only if the 
relative orientation is appropriate, otherwise they have to 
move apart. Since in a frustrated loop the spins cannot 
satisfy all interactions, in this model particle configura- 
tions in which a frustrated loop is fully occupied are not 
allowed. The frustrated loops in the model are the same 
of the spin glass model and correspond in the liquid to 
those loops which, due to geometrical hindrance, cannot 
be fully occupied by the particles. 

Another possible interpretation of the model is that 
m is the occupation variable of the i-th cell and Si in- 
dicates the position of the center of mass of the particle 
inside the cell. Si can in principle assume many differ- 
ent values corresponding to the coordinates of the center 
of mass inside the cell. The interaction should take into 
account that not all pair of internal degree of freedom in 
two neighbors cells are allowed. In the frustrated lattice 
gas model for simplicity we drastically reduce the inter- 
nal degree of freedom to only two values {Si — ±1) and 
introduce the interaction e,j = ±1 between two neighbor 
spin variables to mimic the local disorder of the medium. 
In the first version of the model we assume that the local 
disorder are quenched. 

In the case J = oo the model has a maximum density 
p max ~ 0.68. It has been shown that there exists some 
density p c ~ 0.62, where the system has a transition of 
the type of 3d p-spin model [12] , with a divergence of the 
static non-linear susceptibility 



A. The quenched model 



XSG 



= tfY^( SiniS i n tf 



(11) 



Recently a lattice model, which has mean field prop- 
erties closely related to those of p-spin models, has been 
introduced [5] in connection with the glass transition. 
This model is a diluted spin glass, which, being consti- 
tuted by diffusing particles, is suited to study quantities 
like the diffusion coefficient, or the density autocorrela- 
tion functions, that are usually important in the study 
of liquids. The Hamiltonian of the model is: 



where the average (• • •} is over the Boltzmann measure, 
while the average [• • •] is over the disorder configurations 

{e«}- 

Here we show the results for the relaxation of the self- 
overlap, which is defined as 

q{t) = ^ E Si{t')m{t')Si{t' + t)m{t' + t), (12) 



-/3H = jy^XajSiSj - l)mnj +^Vn,, 



(«) 






(10) 



where Si = ±1 are Ising spins, 



0, 1 are occupation 



variables, and e^ = ±1 arc quenched and disordered in- 
teractions. 

This model reproduces the Ising spin glass in the limit 
p — > oo (all sites occupied, m = 1). In the other limit, 
J — ► oo, the model describes a frustrated lattice gas with 



and for the dynamical susceptibility 

X(t) = N[(q(t) 2 )-{q(t)) 2 } 



(13) 



where the average (• • •) is done on the reference time t' . 
In Fig. 1 it is shown the relaxation functions (q{t)) for 
a system of size 16 3 for densities between p = 0.58 and 
0.62. Each curve is obtained averaging over a time inter- 
val for t' of 6 x 10 6 -8 x 10 7 Monte Carlo steps, and finally 



averaging the results over 16 realizations of the disorder. 
Note that there is no sign of a two step relaxation. The 
long time tail of the functions can be well fitted by a 
stretched exponential form, with an exponent [3 strongly 
dependent on the density, which tends to very low values 
(3 ~ 0.2 at high density. In Fig. 2 it is shown the dy- 
namical susceptibility %(i) for the same size and values 
of density of Fig. 1. Note that %(£) grows monotonically 
and has no maximum at finite time. The asymptotic 
value x(°°) corresponds to the static susceptibility (11), 
and therefore has a divergence at the density p c ~ 0.62. 



B. The annealed model 



We have studied the frustrated lattice gas model (10) 



in the case where the interactions q 



±1 are annealed. 



When we evaluate the partition function of the model, 
we must consider in this case not only the Si and Hi, but 
also the e^ as dynamical variables. Thus, summing over 
the Cij and Si we obtain, apart from an irrelevant factor, 



Z= £ e 

{rii.Si.eij} 



-I3H 



j2 e /3ffetf ; w here 



-(3H eS = -K ^2 n i n j 






(14) 



and K = — In [(1 + e~ 2J )/2] . Therefore the static prop- 
erties of the model are equal to those of a lattice gas 
with a repulsion between nearest neighbor particles, and 
with no correlation between spins, (SiSj) — Sij. With 
the change of variables m = |(1 + <7j), where m = ±1 
are Ising spins, this Hamiltonian can be written as the 
Hamiltonian of an antiferromagnetic Ising model with 
an effective temperature T c g = AK~ X , which is always 
greater than the critical temperature of the 3d antiferro- 
magnetic Ising model T c ~ 4.5. Therefore we can con- 
clude that the model (14), and then also the model (10) 
with annealed interactions, does not present any thermo- 
dynamic transition. In the following we consider always 
the model with J = oo. 

We assume a dynamics for the variables e^ with a ki- 
netic constraint, namely e^ can change its state only if 
the sites i and j, and all their nearest neighbors, are 
empty; in this way the accessible states to a given parti- 
cle may change only if a wide enough region of the sys- 
tem around it rearranges itself. We expect that, as the 
temperature decreases, the disorder due to the local en- 
vironment changes so slowly that the interactions behave 
more and more as frozen playing the role of "self-induced 
quenched" variables. 

In order to generate an equilibrium configuration at 
a given density we simulate the model without any dy- 
namical constraint. In this case we can thermalize the 
system even at high density. Once an equilibrium con- 
figuration is obtained, we consider a diffusive dynamics 
for the particles while the interactions evolve with the 



kinetic constraint, as described before. In conclusion the 
simulations are done in the following way: 

1) one starts from an equilibrium configuration ob- 
tained at some density p; 

2) at each step of dynamics an interaction ey is ran- 
domly chosen and is changed if the sites i and j, 
and all their nearest neighbors, are empty; 

3) a particle (occupied site) on the lattice, one of the 
coordinate directions, and a final state of the spin 
Si are randomly chosen; 

4) one tries to move the particle to the nearest neigh- 
bor site in the chosen direction. The particle can 
move if two conditions are both satisfied. First, the 
destination site must be empty. Second, the spins 
of the particles that are nearest neighbors of the 
destination site, must satisfy the interaction with 
the spin of the chosen particle. If the particle can- 
not move in the chosen site, then the move is re- 
jected; 

5) the clock advances one unit of time. 

During this dynamics we have evaluated relaxation func- 
tions and dynamic non-linear susceptibility Note that as 
density grows, the relaxation time gets longer and longer, 
and gets longer than our observation time (which is be- 
tween 10 7 and 10 8 for a system of size 16 3 ) at a density 
approximately p ~ 0.63. 

In Fig. 3 we show the relaxation functions of the self- 
overlap (12), for a system of size 16 3 , for various densities 
between p = 0.52 and 0.61. Each curve is obtained av- 
eraging over a time interval for t' of 3 x 10 6 -10 8 Monte 
Carlo steps. Observe that for high density the relaxation 
functions clearly develop a two step relaxation, signal- 
ing the existence of two well separated time scales in the 
system. We interpret the first short time decay of the 
relaxation functions as due to the motion of the particles 
in the frozen environment, which on this time scales ap- 
pear as quenched, while the second decay is due to the 
evolution of environment, and final relaxation to equilib- 
rium of the system. The long time tail of the relaxation 
functions is well fitted by a stretched exponential form 
(1), where the exponent (3 depends very weakly from the 
temperature (it is constant within the errors) and ranges 
between (3 = 0.4 and (3 = 0.6. In Fig. 4 we show the time- 
temperature superposition of the relaxation functions of 
the overlap, for densities between p = 0.58 and p = 0.61. 
We tried to fit the intermediate time part, corresponding 
to the plateau, of the relaxation function of the overlap 
for density p = 0.61, with the function predicted by the 
MCT (in a simplified form) 



(q(t)) =f + Af 



Bt b 



(15) 



where the fitting parameters are /, A, B and A, while a 
and b are given by the relation (2). The result is shown 



in Fig. 5, where the full line is the fitting curve with 
a = 0.339 ± 0.002 and b = 0.69 ± 0.01. 

In Fig. 6 we show the dynamical non-linear susceptibil- 
ity (13) for the same size and values of the density of Fig. 
3. It has the same behavior of the p-spin model in mean 
field and of the molecular dynamics simulation of the 
Lennard- Jones binary mixture [3], namely a maximum 
x(t*) that seems to diverge together with the time of the 
maximum t* , when the density grows. For the highest 
density , the maximum of x(t) decreases, possibly due to 
finite size effects, too short observation time, or a change 
in the dynamics above some critical density. This fact is 
observed also in MD simulations of Lennard- Jones liquids 
[4]. We obtain that the maximum x(**) as a function of 
the density can be fitted quite well (taking out the last 
three points, where presumably a rounding of the diver- 
gence takes place) by the power law x(i*) ex (p c — p)~ a , 
with p c = 0.66 ± 0.01 and a = 3.6 ± 0.2. At very long 
times x(t) decays to the equilibrium value, which is sim- 
ply x(oo) = p 2 . 

To make a more direct comparison with MCT, we have 
evaluated, on a cubic lattice of size 8 3 , the autocorrela- 
tion functions of the density fluctuations 
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(6 Pq (t' + t)6 Pq (t')) 

(\S Pq \ 2 ) 



(16) 



where the average (• • •} is performed on time t', the den- 
sity fluctuation of wave number q is defined by Sp q — 
p q - (p q ), and 



p q {t) 



^ e -i,.r,(t) ) 



(17) 



where Ti(t) is the position of the i-th particle at time t 
and n is the particle number. 

Because of periodic conditions the allowable values of 
q on a cubic lattice are of the form 



2tt 
Q = —{n x ,n y ,n z ), 



L 



(18) 



where 



= 1. 



n x ,n y ,n z = 1 • • ■ ■§ are integer values. 

In Fig. 7, 8, 9 we show the results obtained at den- 
sities between p — 0.380 and p = 0.602 respectively for 
q = (tt/4,0,0), (tt/2,0,0) and (tt,0,0). Each curve is 
obtained averaging over a time interval for t' of 10 6 - 
10 7 steps and finally averaging the results obtained by 
32 — 128 different simulations. As we can see in figures, 
at low density the autocorrelation functions relax to zero 
with a one step decay; on the other hand, as the density 
increases we can recognize the two-step decay character- 
istic of glass forming systems. 

As we have said in Sect. II, if Q q (t) satisfies the pre- 
diction of the MCT, R q (t) (7) must be independent on q 
in the (3- regime. We have evaluated R q (t) (with t" ~ 400 
and t' ~ 1.6 • 10 5 ) at some densities for all values of q 
considered here and we have obtained that R q (t) is inde- 
pendent on q on a large time interval (see Fig. 10). In 



agreement with this result we find that, after the initial 
transient, the correlators <& q (t) are well fitted by a power 
law: 



; " + ""l£ 



(19) 



and, in the intermediate time region, by a von Schwcidlcr 
law: 



fq - K 



(20) 



where the exponents a and b are independent on q. 

In Fig. 11 and 12 we show [3> 9 (£) — f q ] as function 
respectively of t and t/r for q = (ir/4, 0, 0) at different 
values of densities. As we can see in figures the curves 
scale for all values of density considered here and the data 
are in good agreement respectively with the power law 
(t/t ) a (with a = 0.66 ± 0.11 and t = 3.6 ± 3.0) and the 
von Schweidler law -{t/r) b (with b = 0.80 ± 0.13). The 
relaxation time t, obtained as fitting parameter from Eq. 
(20), is an increasing function of density, well represented 
by a power law (p c — p)^ 1 , with p c = 0.66 ± 0.03 and 
7 = 1.5 ± 0.2 (see Fig. 13). However the values of the 
exponents a and b obtained in this case do not satisfy the 
relation (2) and are not in agreement with the exponents 
obtained for (q(t)). 

In Fig. 14 we show the correlators & q (t) as functions 
of the rescaled times t = t/r a ; as we can see, the curves 
coincide at large t, for all values of density considered 
here, with a common master curve $ g , well fitted by 
a stretched exponential function (for q = ("7r/4, 0, 0) we 
obtain (3 ~ 0.45). This result is consistent with the pre- 
dictions of the MCT concerning the a— relaxation. 

Finally as the density grows the equilibrium system 
smoothly evolves towards an "ordered" state, analogous 
to the crystal state (the system at high enough density 
can reach the equilibrium state only if frustration is re- 
duced). As a consequence of this fact, the similarity be- 
tween the annealed model and a supercooled glass form- 
ing liquid fails at high density. 



V. CONCLUSIONS 

The frustrated lattice gas model in the quenched ver- 
sion presents a thermodynamic transition at a critical 
density, where the static non-linear susceptibility di- 
verges. The annealed model, which does not present any 
thermodynamic transition, consequently does not show 
any critical behavior of the static susceptibility; on the 
other hand we observe an "apparent" divergence of the 
dynamical susceptibility at a value of density p c (where 
the structural relaxation time also diverges). We sug- 
gest that this behavior is due to the fact that at short 
enough time the disorder can be considered as quenched. 
Moreover the similarity between the annealed model and 



the supercooled glass forming liquids suggests that also in 
these systems a similar mechanism may be responsible for 
the critical behavior of the dynamical susceptibility, and 
one might speculate that if one could in special systems 
freeze some degree of freedom one could find a behavior 
similar to systems with quenched disorder. 
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FIG 2. Dynamical susceptibility in the quenched model, 
for the same system size and densities of Fig. 1. 
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FIG. 3. Relaxation functions of the self-overlap in the an- 
nealed model, for a system of size 16 3 and densities p = 0.52, 
0.53, 0.54, 0.55, 0.56, 0.57, 0.58, 0.59, 0.60, 0.61. 
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FIG. 4. Time temperature superposition principle for the 
relaxation functions of the self-overlap, for densities p = 0.58, 
0.59, 0.60, 0.61. The fitting function is a stretched exponential 
with exponent f3 = 0.5. 




FIG. 7. Correlation functions of the density fluctuations, 
$ 9 (i), for q = (-7r/4, 0, 0) at densities (from bottom to top) 
p = 0.380, 0.490, 0.543, 0.584, 0.602. 
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FIG. 5. Fit of the intermediate time part of the relaxation 
function of the self-overlap, for density p — 0.61, with the 
fitting function / + At~ a — Bt b , where the fitting parameters 
are /, A, B and A, and a and b are given by the relation (2). 



FIG. 8. Correlation functions of the density fluctuations 
$ q (i) for q — (7r/2, 0, 0) at densities (from bottom to top) 
p = 0.380, 0.490, 0.543, 0.584, 0.602. 
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FIG. 6. Dynamical susceptibility in the annealed model, 
for the same system size and densities of Fig. 3. 
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FIG. 9. Correlation functions of the density fluctuations 
$ g (i) for q = (7r,0,0) at densities (from bottom to top) 
p = 0.380, 0.490, 0.543, 0.584, 0.602. 
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FIG. 10. R q (t) for q = (tt/4,0,0), (tt/2,0,0) and (tt,0,0) 
at density p = 0.584. 



FIG. 13. The relaxation time r as function of density p. 
The full line is the power law (0.66 — p)~ . 
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FIG. 11. \<b q - f q ] as function of time t for q = (tt/4, 0, 0) 
at p = 0.548, 0.573, 0.602, 0.615. The full line is the power 
law (t/3.6)^ ' 66 . 



FIG. 14. $ g as function of t/r a for q = (tt/4,0,0) at 
p = 0.467, 0.490, 0.521, 0.543, 0.573, 0.584. The full line 



is the stretched exponential function e 



-(i/Vc) - 45 




FIG. 12. [$, - /,] as function of t/r for q = (tt/4, 0, 0) at 
p = 0.548, 0.573, 0.602, 0.615, 0.625. The full line is the von 
Schweidler law -(i/r) ' 80 . 



